Demonstration of quantum Zeno effect in a superconducting phase qubit 
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Quantum Zeno effect is a significant tool in quantum manipulating and computing. We propose its observation 
in superconducting phase qubit with two experimentally feasible measurement schemes. The conventional 
measurement method is used to achieve the proposed pulse and continuous readout of the qubit state, which are 
analyzed by projection assumption and Monte Carlo wave-function simulation, respectively. Our scheme gives 
a direct implementation of quantum Zeno effect in a superconducting phase qubit. 
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Quantum Zeno effect (QZE), proposed by Misra and Sudar- 
shan in 1977 ll|], is a paradigm showing that quantum physics 
is counter-intuitive. It predict that if the state of a unstable 
or oscillating quantum system is measured frequently to see 
whether it still stay at a initial state, transitions from the initial 
state to other states will be suppressed or even inhibited. Since 
then, many exciting progresses have been made both theoret- 
ically and experimentally. In the theoretical side, physicist 
interpret it with wave-function collapse assumption in early 
days J3, which is shown to be not necessary 0]. Later, 
it was generalized in a few different ways. Concerning the 
measurement, it can not only retard incoherent decay but also 
coherent Rabi oscillation. On the other hand, for unstable sys- 
tem, frequent measurements may even also enhance the decay 
rate under some conditions, which is the so-called quantum 
anti-Zeno effect H-lzI]- As to the readout aspect, one can adopt 
pulse or continuous measurements JsllHl- In the experimental 
side, QZE have been demonstrated in many systems, such as 
trapped ions L3], optical lattice Bose-Einstein condensate 
01, microwave cavity lHoh . etc. 

Studying QZE is very important. Beyond the interest 
of fundamental physics, it has many practical applications. 
These includes reducing decoherence in quantum computing 
fTTI - O . efficient preservation of spin polarized gases [14], 
keeping system stay in object subspace W15I1 . There are inter- 
esting explorations of applications in superconducting qubit 
systems, e.g., generation of entangled state MM and imple- 
mentation of quantum switch W17I1 . Recently, the possibil- 
ity of observing QZE in superconducting qubits is proposed 
0j||2(J. However, demonstration of QZE in superconducting 
system is very difficult because of the lacking of competent 
measurement method. Conventionally, it was observed with 
quantum non-demolition (QND) readout. In circuit Quan- 
tum Electrodynamics (CQED), state of the qubit could be im- 
printed on the cavity field state in a QND readout, but the 
signal-to-noise ratio is so low that we must repeat consider- 
able times to complete the QND readout. Thus, a recent exper- 
iment [21] demonstrate QZE qualitatively in CQED can not 
guarantee its practical applications due to the noise. Here, we 
suggest experimental feasible schemes to demonstrate QZE in 
a superconducting phase qubit with both pulse and continuous 
measurement strategy instead of QND readout. 

Since the quadratic decay behavior (prerequisite for QZE) 
in the initial decay stage of the qubit excited state has not 
been observed yet, suppressing of energy relaxation in super- 





FIG. 1 : Energy levels in the potential well when measurement pulses 
is off (a) and on (b); c, indication of the bias of the phase qubit, 
including microwave driving I MW (upper) and dc bias current with 
measurement pulses (down). 



conducting qubit is not accessible technically. Thus, we here 
focus on another case of QZE, i.e., suppressing the unitary 
evolution of the phase qubit. There are at least two feath- 
ers that differentiate our schemes from those implemented in 
other systems. Firstly, the measurement method used here 
is the so called selective measurement instead of QND read- 
out which is still a big challenge to realize continuous mea- 
surement in superconducting qubits system. Secondly, our 
schemes are immune from the relaxation of the qubit by us- 
ing an appropriate initial state. That is rather necessary in the 
context of the very short energy relaxation time of supercon- 
ducting qubits. It should be noticed that, besides the function 
of demonstrating the basic phenomenon of quantum mechan- 
ics, our proposal can lay a foundation for the applications of 
QZE in quantum information processing, e.g., Ref. [UrSll . 

Superconducting phase qubit usually consists of a large 
current-biased Josephson junction (JJ). When the bias cur- 
rent approaches its critical current, there exist several no- 
degenerate energy levels in each well of the washboard po- 
tential of the qubit, see Fig. 1(a). The lowest two states act 
as a qubit, which is the so-called phase qubit. Experimentally, 
it can be easily controlled by bias cunrrent containing both dc 
and microwave components [ 22h . 

An advantage of phase qubit over other types of supercon- 
ducting qubit is its built-in readout. It relies on the possibility 
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that the qubit states in the potential can tunnel through the 
potential barrier into continuum outside. The tunneling rate 
of one level usually differ dramatically from the other one at 
least two orders of magnitude fl23ll . So the ground and excited 
state can be mapped to no-tunneling and tunneling case, re- 
spectively. Experimentally, one can lower the barrier so that 
the excited state of the qubit can tunnel through the barrier 
quickly but the ground state can't, see Fig. 1(b). Therefore, 
one can add a pulse to the bias 112411 . see Fig. 1(c), so that the 
height of the barrier only fitting the excited qubit state. 

Now, we begin our pulse measurements scheme. When op- 
erating the phase qubit for quantum gates, the bias is tuned to 
an appropriate value so that there are three or four levels in 
the potential well as shown in Fig. 1(a). In this case, neither 
of the qubit states could tunnel outside. Drive the phase qubit 
with a resonant microwave, the qubit will oscillate between 
the ground and excited state with a period of 2n/£l. If the 
initial state is the ground state, after half a period, the qubit 
is driven to the excited state. To demonstrate QZE, we su- 
perpose a series of short uniform measurement pulses to the 
qubit bias, see Figl. (c). Because the pulse during time r 
is much smaller than the oscillating period, we consider that 
the probe pulse is instantaneous. To be more specifically, in 
half a period of Rabi oscillation T = n/Cl, there are n evenly 
distributed pulses with the time interval as St = T/n. 

According to the Hamiltonian, we can calculate straightfor- 
wardly the population of the ground state |0) at t = St, before 
the first pulse, as P^(St) = cos 2 (n/2n) . After the first probe, 
the probability of no-tunneling Pj, i.e, the probability of the 
qubit collapsing to ground state, equals to P^{St). So, after all 
the ni\\ probes, the survival probability of the initial ground 
state |0> is P n = (pj)" = [cos 2 in/In)]" . When n » 1, mak- 
ing the proximation of cos(^) « 1 - (nlln) 2 and using the 
relation lim(l - x/n) n = e~ x , one gets 
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We have plotted the survival probability with both expressions 
in Fig. 2, which shows that the approximation to exponen- 
tial function is perfect when n > 10. Obviously, with the 
increasing of n, the survival possibility of initial state tends to 
1, which is one kind of QZE. 
Although it is similar with the experiment with trapped ions 
there has significant difference between them. In their 
experiment, they use a series of QND measurements during 
the Rabi oscillation and measure the final state of the trapped 
ions at the end of a half oscillating period. Therefore, even 
if some measurement results are not the initial state, the qubit 
still evolve according to Hamiltonian after the measurements. 
Thus, there exist a small probability to return to the initial state 
at last. Instead, we employ a selective measurement approach 
to obtain the probability that all the probes get the same result, 
i.e, the initial ground state. This is achieved by the fact that 
if the result is other than the initial state, the JJ will switch to 
a non-zero voltage state, which means that the state of qubit 
will be destructed and stop evolve after the probe. It should be 




FIG. 2: Survival probability of |0) in a half period of Rabi oscillation 
vs. number of probes. Red dots is plotted with PJ = [cos 2 (7r/2«)j , 
and solid line with P n Q = exp(-7r 2 1 An). The approximation to expo- 
nential function is perfect when /i > 10. 
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FIG. 3: Survival probability of |0> vs. time. Blue solid line is 
showing Rabi oscillation without any decay. Green and black line 
is plotted according to cos 2n (£^/2) with St = n /50£2, n/ 100Q, re- 
spectively. Purple circle and red dot denote exp((-Q 2 St/A)t) with 
corresponding St as above. 



noted that our scheme exhibit what Misra and Sudarshan first 
called QZE [H. 

It is well known that when the interval time among the mea- 
surement pulses is very small, the survival probability of the 
initial state reduces exponentially with the increasing of time. 
Similarly, we can easily get 



Po(t) = Po * exp 



QrSt 



(2) 



with t = nSt. We have plotted P (f) with St = tt/50Q, tt/ 100ft 
in Fig. 3. Instead of normal Rabi-type oscillation, the initial 
state |0> decay exponentially with an effective characteristic 
time t c 
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It is a very important characteristic parameter showing that 
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to what extent the probes suppress the state transition. When 
St = 7r/50Q, t c = 200/7rQ, which is much larger than the 
characteristic time of Rabi oscillation 1 /Q, we can also see 
from Eq. © that the decay time is inverse function of St. 

Then, we move to the feasibility of the above pulse mea- 
surements scheme. Theoretically, the more frequently a qubit 
state is observed, the more likely it will be inhibit to transi- 
tion to other state from the initial state. However, in practice, 
there are three main obstacles stopping us from beating the 
QZE limit. Firstly, measurement fidelity is always lower than 
unit; secondly, each measurement is not instantaneous but in- 
evitably lasts for a finite period of time; and finally, the finite 
decoherence time of the qubit. 

The measurement of phase qubit state is achieved by using 
its macroscopic quantum tunnel. The imperfect fidelity is in- 
duced from the finite ratio of the tunneling rates of |1> and |Q) 
states. It's believed that the ratio is typically around 200 Il23ll . 
During the measurement pulse, the tunneling possibility from 
the excited state is close but a little bit lower than 1 , while that 
of the ground state is a small but nonzero quantity. However, 
this measurement have single shot readout fidelity up to 96% 
theoretically 112311 . which is the highest among all the known 
readout approaches of superconducting qubits. 

The other bothering factor is that the readout can't be ac- 
complished instantaneously. In a measurement process, the 
added probe pulse alters the level structure of the qubit, mak- 
ing it detuning from the driving microwave. The question how 
to judge the suppression of the oscillation comes from QZE or 
from the reduce of the effective driving time is unavoidable to 
any experimental scheme of demonstrating QZE. Actually, for 
larger measuremnt times n, the sum of the measurement pe- 
riods is not negligible compared to the duration of the whole 
process. Therefore, part of the decrease in the transition prob- 
ability is due to the decrease in the time during which the qubit 
is resonant driven. For an extreme case of n = 100, the sum 
of the measurement periods is 50% of the total time T. Even 
for this case, the survival probability of the initial state is as 
high as 97%, which is much higher than that of sole resonant 
driven. So we can safely conclude that the suppression of the 
oscillation mainly comes from QZE. 

Finally, the decoherence including relaxation and pure de- 
phasing of the qubit is usually considered as bottle-neck for il- 
lustrating QZE. From now on, we would clarify why it can be 
neglected in our scheme. On the one hand, pure dephasing can 
be ignored after noticing that the interval between two nearest 
measurements is very short compared with the dephasing time 
which could be as long as hundreds of ns. On the other hand, 
the lifetime of an average phase qubit is T\ = 600 ns. If T 
can be much smaller than T\ , to avoid the decay of the excited 
state, then the experiment will be able to carried out easily. 
But that is not necessary for observing QZE. The reason is 
that before each probe, the average population of excited state 
is much lower than 1, i.e., 

1 r T 

Pi = - sin 2 (Q;/2)<^ = QV/12 « 1. (4) 

T JO 

To be more specifically, P\ 0.07 for n = 32, i.e., each probe 
will project the qubit sate to the ground sate subspace with 



high probability. Therefore, the qubit state can only has a 
small probability to excite to the excited sate, which means we 
will have a much longer effective lifetime for the excited state 
in our scheme. If we conservatively choose T = T\ = 600 ns 
and each probe pulse lasts t p = 3 ns, then the ratio between 
the measurement time and the total oscillating time is T/t p = 
200. In the trapped ions experiment, the quantity is about 100. 
With a larger ratio, one can implement more probes within a 
half period of Rabi oscillation. So, we strongly believe that 
QZE can be verified definitely in phase qubit with our pulse 
measurement method. 

Next, we propose to demonstrate QZE in a phase qubit by 
continuous measurement. Theoretically, Heisenberg uncer- 
tainty principal limits how frequently a measurement can be 
performed. However, one could also adopt continuous mea- 
surement approach to observe QZE. The main difference of 
the continuous measurement approach from the above pulse 
measurement scheme is that the bias is fixed to only allow 
the excited state to tunnel outside during the oscillation of the 
qubit. This is reasonable since the tunneling rate of the excite 
state is two orders of magnitude larger than that of the ground 
state i23ll . Furthermore, the tunneling rate of the ground state 
is also smaller than the Rabi frequency of the qubit in our pa- 
rameter figuration. Therefore, we do not need to take the rare 
tunneling event of the ground state into account. The system is 
initially prepared in the ground state; a resonate microwave is 
driving the phase qubit between the ground and excited state. 
The life time of the excited state due to spontaneous decay is 
much longer because of quantum tunneling is very fast. Ne- 
glecting the spontaneous decay term, we get a Hamiltonian 
describing this dissipative system in interaction picture is 



Before going into QZE, we would like to discuss how con- 
tinue measurement works. In our case, the excited state has 
a tunneling rate T, but the ground state can't tunnel, which 
means in a short time interval St, the excited state tunnel with 
the possibility of TSt. If a tunneling count, we know the qubit 
state before tunneling is the excited state; otherwise we can't 
discern the ground and excite state, but what we can get from 
the interrogation is that the qubit is more likely in ground state 
at the end of St than at the beginning. This point is the the es- 
sential of Monte Carlo wave-function method, which is devel- 
oped for simulating open system Il25l - l27ll . Below we use this 
method to show QZE and compare it to the analytical result. 

Back to the Hamiltonian in Eq. ([5]). If T = 0, the qubit is 
absolutely populates |1) at t = T = n/2£l. Now, we suppose 
Q <sc T. The non-Hermitian Hamiltonian can be simulated 
with Monte Carlo wave-function method. The procedure can 
be summarized as follows: (1). Discretize the time interval T 
by a very small time step St. (2). Determine the probability 
of tunneling P = TSt\(l\i//)\ 2 , choosing St to make sure P <?c 1. 
(3). Obtain a random number r distributed uniformly between 
zero and one, and compare it with P. (4). If r < P, there 
is a tunneling, the system switch to finite voltage state, and 
this run is end. Then start next run from step 1. If r > p,no 
tunneling takes place, the qubit evolves under the influence of 
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FIG. 4: No-tunneling counts vs tunneling rate of the excited state. 
Monte Carlo wave-function simulation parameters are Q = 2n x 1 
MHz, T = tt/2Q and T e [50, 500] MHz. The blue circles are data 
from simulation, and the red solid line is plotted according to the 
analytical expression in Eq. (|8]). 



the non-Hermitian Hamiltonian described by Eq.© and the 
qubit state at the end of St is 

iHSt iHSt 
m + St )) = (1 - )|^(0)/||(1 - —M(t))\\ (6) 
n n 

where we have approximately expand the the evolution oper- 
ator to first order of St. (5). Repeating the process, we can get 
a trajectory of the qubit state. 

It is obvious that if no tunneling appears in the whole period 
T, the qubit state will totally stay in the state of |0> at the end 
of the probed oscillation. So, the survival probability of a ini- 
tial state is same as that of no tunneling during the period. We 
simulate many times for a certain T to obtain survival proba- 
bility Po. Furthermore, we have studied the relation between 
Po and T as shown in Fig. 4. 

The Hamiltonian in Eq. ([5]) can also be solved analytically. 
The evolution operator of the dissipative two-level system has 
the form of 



U = e~ iHlt e- 



h • a 

cosh(ht) - i — - — sinh(/z0 



where h = V074) 2 - Q 2 and we have assumed that T/4 > Q. 
If the initial state is |0>, then the survival amplitude is function 
of time and has the form of 



A (0 



<0|£?-^|0> 



= e 



cosh(/z0 + — sinh(/z0 
2h 



(8) 



(7) 



We have also plotted the survival probability of the initial state 
with this analytical expression in Fig. 4. We can see the 
Monte Carlo wave-function simulation is agree with the an- 
alytical result perfectly. More importantly, they both imply 
QZE as explained in the following. With the increasing of the 
excited state tunneling rate, the no-tunneling counts approach- 
ing to 1000, which is the total simulation runs' number. We 
conclude that the survival probability of the initial state |0) 
is more strengthened with larger tunneling rate of the excited 
state. That is the essential meaning of QZE of a system mea- 
sured continuously. 

Additionally, we can also investigate the relation of pulse 
and continuous scheme in the future experiment. Theoreti- 
cally, it could be proved that if the effective decay times of 
the initial state in the two schemes are the same, the interval 
between two sequential measurements in pulse scheme and 
decay rate of the excited state in continuous scheme should 
satisfy the relation [8D of St • T = 4. That is an important 
relationship between the two schemes. 

We have proposed two schemes to observe QZE in a su- 
perconducting phase qubit: pulse and continuous measure- 
ment schemes. They are easy and feasible for up-to-date tech- 
nique. Our result show that QZE can be demonstrated with 
the schemes more clearly than the trapped ions experiment. 
Since QZE is essential in the implementation of quantum in- 
formation, the generalization of the proposed scheme to other 
kinds of superconducting qubits |28, 29] is desirable. 
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